Numerical simulations of a ballistic spin interferometer with Rashba spin-orbital 

interaction 
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We numerically investigate the transport behavior of a quasi one-dimension (ID) square loop 
device containing the Rashba spin-orbital interaction in the presence of a magnetic flux. The 
conductance versus the magnetic field shows the Al'tshuler-Aronov-Spivak (AAS) and Aharonov- 
Bohm (AB) oscillations. We focus on the oscillatory amplitudes, and find that both of them are 
strongly dependent on the spin precession angle (i.e. the strength of the spin-orbit interaction) and 
exhibit no-periodic oscillations, which are well in agreement with a recent experiment by Koga et 
al. @- However, our numerical results for the ideal ID square loop device for the node positions 
of the amplitudes of the AB and AAS oscillations are found to be of some discrepancies comparing 
with quasi- ID square loop with a finite width. In the presence of disorder and taking the disorder 
ensemble average, the AB oscillation in the conductance will disappear, while the time-reversal 
symmetric AAS oscillation still remains. Furthermore, the node positions of the AAS oscillatory 
amplitude remains the same. 

PACS numbers: 71.70.Ej, 72.20.-i, 73.23.Ad 



The control over the transport properties of electron 
spins has gained much attention recently 0,0 01- A new 
sub-discipline of condensed matter physics, spintronics, 
is emerging rapidly and generating great interests. The 
promising application of spintronics widely lies in, e.g. 
information technology the spin electron apparatus [jj, 
and so on. Utilizing the spin orbital (SO) interaction 
to manipulate the spin degrees of freedoms of electrons 
has been advanced. Several spin devices have also been 
theoretically designed. Experimentally, the strength of 
Rashba SO interaction has successfully been controlled 
by the applied gate voltage or with some specific design 
of heterostructures 0, . 

Very recently, Koga et al. used a nanolithographically 
defined square loop array in the two dimensional electron 
gases to experimentally demonstrate a gate-controlled 
electron spin interference 0, |j| . They observed that the 
amplitude of the Al'tshuler-Aronov-Spivak (AAS) type 
oscillation 01 of the conductance versus the magnetic 
field B depends strongly on the gate voltage, or the sheet 
carrier density. This implies that electron spin preces- 
sion angle 6, caused by the Rashba SO interaction 
is gate-controllable, and can be tuned more than 0.757T. 
In fact, to experimentally realize a large controllability 
of the spin precession angle is a very important issue in 
spintronics, because it is indispensable for realizing some 
spin devices, e.g. the spin field effect transistor y|- 

Before the experimental work ||, some theoretical 
works have investigated the transport behavior of a ring 
with the SO interaction [Til Hii . Il4| . In particular, a 
theoretical prediction of the amplitude of the AAS oscil- 
lation has been made by the same group 15J. Although 
both the experimental and theoretical results are similar, 



there are some discrepancies. The theoretical system is 
an ideal perfect one dimension (ID) square loop with only 
one terminal contacted outer, for which the backscatter- 
ing probability was calculated. In the experiment, the 
device is a quasi-lD square loop (i.e. each arm of the loop 
has a finite width) with multi terminals, and the trans- 
port conductance instead of the backscattering probabil- 
ity was measured. The authors argued that the results 
for both systems should be similar, however, no through 
investigation was given. Moreover, in the experimental 
device, the first order oscillation is the Aharonov-Bohm 
(AB) oscillation with the period tfio (cf>o = h/e), which is 
never studied. 

In this paper, we numerically simulate the experimen- 
tal set-up by investigating the amplitude of the AB and 
the AAS oscillation. We consider a quasi-lD square loop 
having the Rashba SO interaction (see Fig. la), which 
represents one of the cells in the experimental setup by 
Koga et al. . The conductance is first calculated from 
the Landauer-Biittiker formula. Afterwards, by apply- 
ing the same data analysis process as in the experiment, 
i.e. the fast Fourier transform (FFT) and the inverse 
fast Fourier transform (IFFT), the amplitude of the AAS 
and AB oscillations of conductance are obtained. The 
numerically calculated node positions of amplitude of 
AAS oscillation in quasi- ID device are of some depar- 
tures (~ O.I71") from the theoretical results for the ideal 
ID device. The difference is larger for a loop with a wider 
arm. Moreover, the node positions of the amplitude of 
the AB oscillation are also studied. Finally, we discuss 
how the node positions affected by the Dresselhaus SO in- 
teraction and the disorder, and find that the Dresselhaus 
SO interaction can strongly shift the node positions, but 
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the disorder has almost no effect on the node positions. 




FIG. 1: Sketchy illustration for a quasi-lD (a) and an ideal 
ID (b) square loop devices. In (a), the potential in the center 
region is assumed to be very high so that there are no electrons 
there. 

Before investigating the quasi-lD loop system, we first 
study an ideal ID square loop model attached with two 
ID leads on its two opposite corners (shown in Fig. lb). 
The goal is to find analytically the node positions of the 
amplitude of AB or AAS oscillation in the two-terminal 
system. The Rashba SO interaction exists only in the 
loop, but absent in the two leads. [lfij An incident wave 
if>i is splitted at the lower left corner. And the counter- 
clockwise (CCW) and clockwise (CW) wave ballistically 
travelling along two sides of square loop merge at the 
right upper corner. ip a describes the output wave func- 
tion of the electron at the right upper corner. If only 
the first order tunneling process is considered, ip a can be 
obtained as: 

4>o = \ (e-^ 2 R^ x (6)R y (0) +e i */ 2 R y (6)R- x (d)) ^.(1) 

where <f> — 2?r$/0o with the magnetic flux <E>. The rota- 
tion operator Rf(9) is defined as:[l5| 



R f (9) = Icos(0/2) - if • o-sin(0/2), 



(2) 



and 9 — 2m*anL/h 2 is the spin precession angle. Here 
L is the side length of the ID square loop and a>u is 
the strength of the Rashba SO interaction. To consider 
that the incident electron is spin unpolarized, the output 
probability of the electron is < ip \i>o >, and it's given 
by 
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< 4>o\i>o > = -+A(9)cos(f>, 



(3) 



where A(9) — l/4(sin 2 + 2cos0) is the amplitude of 
the AB oscillation. Then the node positions (marked 



by 9*) can be obtained easily, and they are 0.6367T and 
1.364-7T, etc. The amplitudes of the higher order oscil- 
lation (including the AAS oscillation) are zero here, be- 
cause the higher order tunnelling process has been ne- 
glected in Eq.(l). 

Next, we study the model of a quasi-lD square loop 
sketched in Fig. la. With the model, no electron exists in 
the center region of the loop, which can be experimen- 
tally realized by the etching technology or the deposited 
metal split gate. Four leads symbolized as hatched re- 
gions in Fig. la are attached to the four sides of the sys- 
tem. Rashba SO interaction and magnetic field only exist 
in the quasi loop, and four leads are ideal without SO in- 
teraction [J^. The width of the lead is N, the channel 
width is M, and the side length for the quasi-lD square 
loop is L=M+N. This quasi-lD model is identical with 
each cell of the experimental setup in the Ref. ||. In 
their experiment, they applied a square loop array to de- 
termine the sheet conductance which is necessary to di- 
minish the AB oscillation and the universal conductance 
fluctuations. 

The Hamiltonian for the quasi- ID model is described 
as a discrete lattice tight binding model in our numerical 
calculations. We choose a symmetric gauge and the vec- 
tor potential A(x,y) = B(—y/2,x/2,0), with the lattice 
spacing as the length unit. Let a,j a (a] CT ) be the annihi- 
lation (creation) operator of an electron with its spin a 
at the lattice site j, and a] = (aL,at,). Then the tight 
binding model Hamiltonian for the system can be written 
as: 
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+ y^^H^jx a \ cr y a j+x - ^jyOjPxaj+y) + h.c] (4) 
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In the above equation, <J x ( y ) are Pauli matrices, and 
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(5) 
(6) 



where t — 2m *a 2 * s ^ n e nearest neighbor hopping element 
with the lattice space a, A = describes the strength of 
the Rashba SO interaction, and r = ±i, ±y. We did not 
consider the Zeeman splitting energy, because the mag- 
netic field in the experiment is so weak that the Zeeman 
term is negligibly small compared with other energies, 
e.g. t, the Rashba splitting energy, etc. The disorder 
potential at each site is not considered right now. Thus 
we set Wj — in the Hamiltonian as energy zero point. 
The tight binding Hamiltonian for the ideal ID square 
loop system is also given with the same approximations. 

The current I p from the lead p (p = 1, 2, 3, and 4) 
flowing into the square loop can be calculated by the 
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Landauer-Buttiker formulae |l7l Il8| : 

1 

where V q is the voltage applied on the lead q, T pq (E) is 
the transmission probability from the lead q to the lead 
p, and E is the electron Fermi energy. The tempera- 
ture is set to be zero, because the thermal energy ksT is 
much smaller than the other energy scales in the experi- 
ment. The transmission probability T pq is determined by 
T pq = Tr[T p G R T q G A } [3113. Here T p = »[£« - Z A ], 
and G R ( A ^ {E) is the retarded (advanced) Green function 
given by: 

G R (E) = (G A (E)y = l — , (8) 

E - H c - £ 

where E^" 4 -* is the retarded (advanced) self-energy of 
lead p and H c is the single particle Hamiltonian given by 
Eq. (4) for the isolated finite-size system at the center. 
In the following numerical calculations, the lead's volt- 
ages V p are set as: V\ — V/2 and V3 = —V/2, i.e. a 
longitudinal bias is added in the x direction. The trans- 
verse lead- 2 and lead-4 act as the voltage probes, and 
their voltages V2 and V4 are calculated from the condi- 
tion I2 = I a = 0. Finally, the longitudinal conductance 
a is obtained as: a — I\/V — —I3/V. In the numerical 
calculations, we let t — 1 to be the energy unit. 



at different spin orbital interaction A are shown in Fig. 2. 
These conductance curves clearly exhibit the existence of 
the AB oscillation with the period <fio and the AAS oscil- 
lation with the period cj>o/2. Since each arm of the loop 
has a certain width, these oscillations are not exactly pe- 
riodic. Their oscillatory amplitudes at B = can be 
figured out by using the FFT. From the graph, we also 
find the alternative changing of peak and dip feature in 
a(B) at B=0 with increasing of the Rashba interaction. 
This phenomenon is consistent with the experimental re- 
sults n. 

Then for each curve of the conductance in Fig. 2 with 
a fixed value of the Rashba strength A, we can employ 
the FFT and IFFT techniques to extract out the ampli- 
tudes of the AB and AAS oscillations at zero magnetic 
field. The data analysis is identical with the experimen- 
tal procedure. Since A = ^ and t = 2 m' a ' 2 > ^ ne R- asn ba 
strength A and the spin precession angle 9 is related by 
9 = ^jA. Fig. 3 shows the amplitudes <r(B = 0) of the 
AB and the AAS oscillations at B = versus the spin 
precession angle 9 for different sizes of the device. The 
node positions 9* for both AB and AAS oscillations are 
indicated in each graph. For comparison, the amplitudes 
of the AB and the AAS oscillations of the exact ID square 
loop system (see Fig. lb) are also calculated by using the 
tight binding model and the Landauer-Buttiker formulae, 
in which the higher order processions and the reflecting 
processions have all been included (see Fig. 4). 
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FIG. 2: The longitudinal conductance u(B) versus the mag- 
netic field B for the quasi- ID device at different Rashba in- 
teraction A, with the parameters iV = 12, M = 6, and the 

Fermi energy E = —2.0. The value of A is indicated with 

2 

each curve. B is in the unit of ^j^-. The plotted curves are 
shifted along the y-axis for comparison. 

The longitudinal conductance a of the quasi-lD square 
loop system with N=12 and M=6 versus magnetic field B 



FIG. 3: Numerical results for the amplitude of the AAS (a 
and c) and the AB (b and d) oscillations at B=0 of quasi- ID 
square loop system versus spin precession angle 6. The system 
size is N=12, M=6 for part (a), (b), and N=20, M=10 for part 
(c), (d). The 8 values of node positions are indicated in each 
graph. The electron Fermi energy is E = -2.0. 

For the ideal ID system, the numerical results ex- 
hibit that the node positions 9* of the amplitude of 
the AAS oscillation are 0.420tt, 0.821tt, 1.183tt, 1.584tt etc 
(see Fig. 4a), and 9* for the AB oscillation are 0.6377T 
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ideal ID system, and the errors can be more than 0.l7r. 
For a loop with wider arm, the discrepancy is larger. For 
example, for the arm's width M = 10, the second node 
position 9* of the AB oscillation is at 1.5527T (see Fig. 3d), 
which has an error of 0.1887T with the corresponding value 
1.3647T of the ideal ID device. 

Therefore, we need to emphasize that if the width of 
the loop arm in the device by Koga et al. has been taken 
into account, then the node positions 9* deviate from the 
ones for the ideal ID device. However, their experimen- 
tal conclusion, that the spin precession angle 9 due to 
the Rashba SO interaction is gate-controllable, remains 
valid. Furthermore the tunable range of 9 could be much 
larger than the value stated in their paper since the node 
positions 9* have larger deviations for quasi-lD devices. 



FIG. 4: Numerical results for the amplitude of the AAS (a) 
and the AB (b) oscillations at B=0 in an ideal ID square 
loop system versus spin precession angle 9 with L=30. The 
9 values of node positions are indicated in the figure. The 
electron Fermi energy is E =-1.7. 



and 1.371-7T, etc (see Fig. 4b). These values are calculated 
at electron Fermi energy E = —1.7, which is close to 
the bottom of the band. Because near the band bot- 
tom, the energy momentum dispersion relation is almost 
quadratic, and the spin precession angle 9 is independent 
with E. So the node positions of 9 remain unchanged 
even at different Fermi levels E. These node positions 
are in agreement with the previous theoretical resultsflij 
or the Eq.(3) for the ballistic system, and the errors are 
within 0.017T for all node positions 9*. This means that 
the higher order tunneling processions (e.g. electron goes 
around the loop multi-times) and the reflecting proces- 
sions have limited effect on the node positions 9*. Fur- 
thermore, our numerical calculations also find that so 
long as the lengths of CW and CCW paths are the same, 
the conductance is identical no matter where the leads 
are connected to the perimeter of the loop . 

Next, let us discuss the numerical results of the quasi- 
lD system (shown in Fig. 3). The amplitudes Aa(B = 0) 
of the AB and AAS oscillations are oscillatory functions 
of the spin precession angle 9, which are similar with 
the ideal ID device. Moreover, it also clearly shows that 
A<j(B = 0) is not a periodic function. This is contrast to 
an ideal ID device. In that case, Aa(B = 0) versus 9 is 
approximately a periodic function with the period tt for 
the AAS amplitude and 2ir for the AB amplitude, like 
A(9) in previous theoretical calculation [l5| and Eq.(3) 
in this paper. But this fact is consistent with the exper- 
imental results [8j. 

Following, we focus on the node positions 9* of the 
oscillatory amplitudes, which have been given in Fig. 3. 
These values of 9* have quite large difference with an 
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FIG. 5: Numerical results for the amplitude of the AAS (a) 
and the AB (b) oscillations at B=0 of an ideal ID square loop 
system versus Rashba dependent spin precession angle 9 with 
L=30. r is the ratio of Dresselhaus versus Rashba interaction 
strength with r = an /otR. The Fermi energy is E=-1.7. 

Up to now, we only consider the existence of the 
Rashba SO interaction in the loop, and the Dresselhaus 
SO interaction has been neglected. If there exists both 
Rashba and Dresselhaus SO interactions, how are the 
node positions 9* affected? In fact, the ID loop having 
two kinds of SO interactions has been investi gat ed by a 
very recent theoretical work by Ramaglia et aL[21| . Their 
results show that Dresselhaus SO interaction strongly 
shifts the nodes of Rashba-dcpcndcnt transmission in an 
ideal ID model. Here we further investigate the ideal 
ID system with these two kinds of SO interaction. The 
tight binding Hamiltonian of Dresselhaus SO interaction 
is given as: 

H D = j v a]cr y aj+$ - I3 ]x a]a x a j+i ) + h.c] (9) 

3 

where (3j T = fl-e l % A " r , [3 = a_o/2a, anda_o is the strength 
of Dresselhaus SO interaction. The total Hamiltonian in 
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the present system is the Hamiltonian in Eq. (4) plus Hr> ■ 
Then the amplitude of AAS and AB oscillation can be 
solved by the same method as in the above. 

In Fig. (5), we present the results of the amplitudes of 
AAS and AB oscillations at B=0 for the exact ID quare 
loop system as a function of Rashba spin precession angle 
9 at different values of Dresselhaus interaction. In this 
graph we could clearly see that with the increasing of 
Dresselhaus SO interaction, the node positions will not 
only vary, even the number of nodes may change. For 
example in Fig. (5a), the number of nodes will decrease 
from 4 to 3 at r — 0.2 (r = od/or), and further to 2 
at r > 0.2. The same situation also happens in the AB 
oscillation amplitude. These results are consistent with 
the analytical results [2l]]. Therefore, if both SO interac- 
tions are presented in the system, using node positions 
to determine the Rashba SO interaction strength is not 
very reliable. However, from the oscillation of the ampli- 
tude AAS or AB versus the gate voltage, it still clearly 
indicates that the strength of SO interaction can be well 
tuned. 
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FIG. 6: Left panel: ensemble averaged conductance as a func- 
tion of magnetic field B at different disorder strength of the 
exact ID square loop system. Right panel: the corresponding 
Fourier transform spectra of the conductance. The location 
of AB frequency and AAS frequency is indicated in the graph. 
The parameters are the Fermi energy E=-1.7, the system size 
L=16, and the Rashba interaction A = 0.1. 



Finally, we study how the node positions 9* are af- 
fected by the ensemble average. Notice that the device 
in the experiment by Koga et al. is a square loop array, 
not a single square loop|8j. So the ensemble average has 
to be investigated in order to make a direct comparison 
with the experimental results. The various ensemble av- 
eraging procedures have been studied in investigating the 
persistent current in a closed mesoscopic ring about ten 
years ago. Those include taking the ensemble average on 
the number of particle, on the chemical potential, or on 



the disorder [22, 23]. The present square loop array is an 
open system, its chemical potential is determined by the 
leads, but the disorder structure for each individual loop 
is random. When giving the disorder structure and the 
chemical potential, the number of electrons in the loop is 
fixed and its value can be non-integer in an open system. 
So in the following we take the average over the disor- 
der ensemble [2^|. To consider the existence of the disor- 
der, the random lattice on site energy Wj in Eq. (4) is 
generated by an uniformly distribution [— w/2, w/2] with 
disorder strength w. All these conductance curves are 
averaged over up to 2000 disorder realizations at w 7^ 0. 
The conductance as a function of magnetic field B at 
different disorder strength w are shown in Fig. (6) left 
panel, with their corresponding FT (Fourier transform) 
spectra on the right panel. With the increasing of disor- 
der strength, both the conductance and the conductance 
oscillation amplitude will decrease. This is due to the 
fact that the electron will be localized at large disorder 
strength. In particular, in the strong disorder case, AB 
oscillation will disappear, but AAS oscillation still exists 
even at the disorder strength w = 1.0, because the AAS 
oscillation is the interference between two time reversal 
pathes and it is independent of the disorder structure. 
These results are in agreement with the experimental re- 
sults and the previous theoretical predictions |8i Il5l] . 

Next, we focus on the node positions 9* after the dis- 
order ensemble average. Fig. 7 shows the ensemble aver- 
aged AAS oscillation amplitude at B=0 of an ideal ID 
loop model at different disorder strength w. w = is not 
shown in the graph, because in that case the oscillation 
amplitude is much larger than w = 0.6 or w — 1.0. It 
is clear that both the node positions of 9 and the oscil- 
lation shape is unchanged versus disorder. These graphs 
demonstrate that with a little strong spin-independent 
disorder and ensemble average, the AB oscillation will 
disappear in the conductance. Meanwhile, the node po- 
sitions 9* of AAS oscillation amplitude still remains the 
same. 

In summary, by using the tight binding model and 
the Landauer-Buttiker formula, we numerically study 
the electron transport through a quasi one dimensional 
square loop with the Rashba spin-orbit interaction. The 
conductance a as a function of the magnetic field B is ob- 
tained, and exhibits the Al'tshuler-Aronov-Spivak (AAS) 
and Aharonov-Bohm (AB) oscillations. These oscillatory 
amplitudes are non-periodic with the Rashba spin pre- 
cession angle 9. These results are in agreement with the 
recent experiment by Koga et al. -8]. To compare with 
an ideal ID square loop device, the node positions of the 
amplitudes of the AB and AAS oscillations have some 
deviations. For a loop with a wider arm, the deviations 
can be quite large. When under the influence of spin- 
independent disorder and ensemble average, the AB os- 
cillation will disappear and only AAS oscillation survives 
in the conductance. In particular, the node positions of 
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FIG. 7: The amplitude of the AAS oscillation at B=0 of an 
ideal ID square loop system versus Rashba dependent spin 
precession angle 6 with L=20. The node positions are identi- 
cal compared with no disorder. The electron Fermi energy is 
E=-1.7. 



the AAS oscillatory amplitude are almost unchanged. 
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